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ABSTRACT 


The  electric  microfield  distribution  in  fluids  of  particles  of 
arbitrary  size,  shape,  charge  and  charge  distribution  is  studied. 

The  Morita- Iglesias  formalism  is  extended  to  include  angular  dependent 
interactions.  The  results  of  the  mean  spherical  approximation  are  discussed 
using  an  exact  relation  for  the  second  moment  of  the  microfield 
distribution,  obtained  in  this  work. In  the  Onsagerian  strong  coupling 
limit  the  sum  rule  for  the  second  moment  simplifies , and  becomes  identical 
to  the  mean  spherical  approximation  result  for  the  second  moment. 


1- INTRODUCTION 


The  recent  revival  of  the  interest  in  the  study  of  the  electric 
microfield  distribution  in  charged  plasmas  has  been  motivated  by  the  use 
of  the  broadening  of  spectral  lines  in  very  hot  plasmas  to  measure 
can  be  used  to  measure  the  temperatures  in  fusion  processes. 

The  study  of  the  fluctuating  electric  microfields  was  pioneered 
by  Holtsmark1  2  who  studied  the  weak  coupling  limit. 

A  considerable  advance  was  the  application  of  the  classical  fluids  methods 
for  evaluating  the  electric  microfield  distribution  in  systems  of  point  charges 
moving  in  a  neutralizing  background,  a  possibility,  that 
was  first  suggested  by  Morita  3  and  formulated  in  terms  of  correlation 
functions  by  Iglesias  *.The  Morita- Iglesias  formulation  (MI  in  what  follows) 
is  a  powerful  tool  for  the  study  of  the  electric  microfield  and  has 
prompted  a  number  of  investigations  on  the  subject5  9 

In  particular,  in  the  last  reference,  Lado9  discusses  the  use  of  the 
Mean  Spherical  Approximation  (MSA)  for  the  evaluation  of  the 
microfield  in  a  one  component  plasma, as  a  starting  point  for  employing 
the  more  accurate  Hypernetted  Chain  (HNC)  integral  equation. 

The  electric  microfield  problem  is  also  of  interest  in 
relation  with  different  dense  systems  other  than  plasmas.  The  electric 
microfield  distribution  in  solutions  in  general,  and  in  ionic  solutions 
in  particular,  is  of  interest  in  various  branches  in  Chemical  Physics,  such  as 
spectroscopy , light  scattering,  liquid  state  chemical  kinetics,  dielectrics  and 
so  on. 

One  of  the  objetives  of  this  paper  is,  precisely  to  study  the  evaluation 
of  the  electric  microfield  distribution  in  systems  composed  of  particles  with 


any  shape,  size  and  charge  distribution. 

For  these  systems  we  consider  the  distribution  of  the  fluctuating 
field  felt  by  a  point  in  an  object  of  finite  size  and  with  an 
arbritrary  multipolar  charge  distribution. 

We  express  it  in  terms  of  the  multipolar  expansion  for  the  interaction 
energy  between  an  imaginary  dipole  placed  at  the  center  of  our  test  particle 
and  the  distribution  of  the  charged  particles  in  the  system. 

Our  result  is  given  in  terms  of  quantities  that  can  be  computed  by 
approximate  liquid  state  theories.  However,  care  must  be  taken  because 
of  the  inconsistencies  which  are  known  to  occur  in  these  theories. 

Because  of  the  difficulties  in  simulating  the  electric  microfield 
distribution  in  a  molecular  fluid.it  is  useful  to  know  exact  relations 
to  estimate  the  accuracy  of  the  approximations.  The  second  moment  of 
the  microfield  distribution  can  be  evaluated  for  a  charged  object  immersed 
in  a  plasma. This  has  led  to  the  very  succesful  APEX  theory  for  plasmas.5 
However, the  second  moment  for  our  general  case  is  much  more  complicated 
as  we  show  in  Section  III,  for  a  hard  charged  object  of  arbitrary 
shape  in  a  fluid  of  charged  hard  spheres  with  point  ions  and  dipoles,  since 
it  involves  three  body  correlations . The  point  is,  that  in  the  Onsager 
limit  10  13,the  three  body  terms  cancel  out,  leaving  a  remarkably 
simple  expression  for  the  second  moment  of  the  electric  microfield. 

In  this  strong  coupling  limit  the  field  of  every  particle  is  completely 
screened  out  by  the  surroundings:  This  limit  is  equivalent  to  Onsager 's 
procedure  of  immersing  the  test  object  into  a  conducting  continuum  media, 
so  that  the  calculation  of  the  internal  energy  becomes  an  exercise  in 
electrostatics . 

We  illustrate  these  points  with  the  application  of  the 


mean  spherical  approximation  to  a  mixture  of  ions  and  dipoles. 

The  layout  of  the  paper  is  as  follows .  In  the  next  section  we 
formulate  the  microfield  problem  and  introduce  our  model. We  also  give 
a  generalization  of  the  MI  expressions  for  the  molecular  case, as  an 
expansion  in  terms  of  electric  multipoles. 

In  Section  III  we  obtain  an  exact  second  moment  relation  for  an 
arbitrary  convex  object  immersed  in  a  mixture  of  hard  spherical  ions  and 
dipoles.  In  Section  IV  we  discuss  the  strong  coupling  limits 
and  make  contact  with  Onsager's  discussion14. 

Section  V  is  devoted  to  mean  spherical  approximation  calculations  on 
ionic  and  dipolar  fluids . Finally , in  Section  VI  we  make  some  general 
remarks  and  conclusions. 


I I -THE  ELECTRIC  MICROFIELD  DISTRIBUTION 


a)  FORMULATION 


Consider  a  mixture  of  s  species  of  particles  and  let  be  the 

number  of  particles  of  species  a.  The  coordinates  of  particle  i  of  species 
a  are 
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(R 
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.  0  > 
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where  R^  is  a  3 -dimensional  vector  which  gives  the  position  of  the  center 

of  particle  i  and  $  ■  ( ) 

a  i 

indicates  the  three  Euler  angles  that  determine  the  orientation  of  the 
particle.  The  electric  field  experienced  at  position  R^  by  an  object 
immersed  in  the  mixture  is 
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Here  X  -  (X  . X  )  denotes  the  configuration  of  all  the 


particles  of  the  mixture. 


The  electric  microfield  distribution  is  in  the  MI  formalism 
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In  these  equations,  the 
average  over  the  configuration 
the  test  object.  We  have  also 
S  function  and  i  1  , is  the 


angular  brackets  denote  the  canonical  ensemble 

& 

of  the  N  +  1  particles  (  N  )  including 

Of  » » 

used  the  integral  representation  for  the  Dirac 
imaginary  unity. 


b)  THE  FOURIER  TRANSFORM  OF  W(e) 


If  the  system  is  isotropic,  A(A)  , 
microfield  distribution,  can  be  rewritten  4 


the  Fourier  transform  of  the 
(A  -  magnitude  of  A  ) 
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where 

r  -  * 

Q(A)-  dX^dX  e  1  (2.5) 

is  the  configuration  partition  function  for  a  system  interacting 
and  the  potential  energy  is  complex 


V(A)-  V  -(i/a)  ?.?0 


(2.6) 


Here  V  is  the  potential  energy  of  the  real  system  of  N  +  1  particles 
particles  when  A  —  0  and  a  ”  (kT)  1  the  usual  Boltzmann  factor. 

Equations  (2.4)-  (2.5)  can  be  expressed  in  the  form  of  the 
free-energy  integral  in  a  coupling-process  type  calculation  3  *  . 


In  A(A)-(l/64 


a 


I 


(2.7) 
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where  the  electrostatic  part  can  be  represented  by 
of  the  form15  16 


a  multipole  expansion 
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Here  R  -|R  | , R  which  are  the  magnitude  and  direction  of  the  vector 

i;  'i 

J>  -4 

Ri^-R( -R^  .  The  orientational  dependence  is  contained  in  the  rotational 
invariants 
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where  the  notation  and  conventions  of  ref.  (17)  for  Wigner's  generalized 
spherical  harmonics  and  3-J  symbols  have  been  used.  The  radial  functions 
in  Eq  (2.11)  are 
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where  £— m+n,and  we  have  omitted  the  greek  subindices  to  alleviate  the 
notation. The  electric  multipoles  are  defined  by 
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q  (r)  is  the  charge  density  distribution  in  the  particle  of  species  a. 

'A 

We  observe  that  the  factor  between  brackets  in  the  integrand  of  Eq . 
(2.7)  can  be  expanded  in  a  similar  way: 
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with 
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The  pair  distributionfunction  has  the  invariant  expansion: 
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The  superscript  (r)  indicates  the  functions  which  involve  the  real 
multipole  distributions  and  (A)  those  of  the  imaginary 

dipole.  Because  of  the  symmetry  of  the  point  dipole  this  last  contains  only 
invariants  with  m  -  1,  n  -  0. 

Substituting  Eqs .  (2.14)-  (2.16)  into  Eq  (2.7)  we  find,  using  the 


orthogonality  properties  of  the  rotational  invariants 
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which  is  the  generalization  of  the  MI  formula  for  a  mixture  of  charged 
particles . 

The  integral  is  proportional  to  the  energy  of  interaction 
between  the  imaginary  dipole  of  the  test  object  and  the  n-pole  of  a 
particle  of  type  a. 

It  should  be  emphasized  that,  the  'origin'  of  the  test  object  (i.e. 
the  reference  point  for  the  multipoles  Eq.(2.13))  is  not  necessarily  at  the 
centerof  the  particle. Our  formalism,  and  in  particular  Eq  (2.18) 
applies  for  the  microfield  distribution  at  any  point  of  a  given  convex  test 
obj  ect 

It  should  be  pointed  out  that  the  test  object  can  be  an  extra  particle 
or  a  particle  belonging  to  one  of  the  species  in  the  mixture. 


Ill -  SECOND  MOMENT  FOR  ION-DIPOLE  MIXTURES 


For  Che  one  component  plasma  of  charged  point  particles  in  a 
neutralizing  background,  Jancovici18  found  that  the  second  moment  of 
the  microfield  distribution  experienced  for  a  test  point  particle  is 

<  <E0.£  >-  kw(p  /0)( z/z9  ) 

(3.1) 

where  e  (e-electron  charge)  is  the  background  charge  density,  z  and  za 
are  the  electrovalence  of  the  plasma  particles  and  of  the  test  particle. 

The  method  used  for  the  plasmas  can  be  extended  to  calculate  the  second 
moment  of  the  microfield  distribution  for  a  hard  convex  object  of  arbitrary 
shape  .charge  z^e, immersed  into  a  mixture  of  arbitrarily  charged  hard  spheres 

O 

We  assume  that  species  a, of  diameter  ,  has  charge  z^e  (  -  Q^(a)) 

-*  1 

and  also  a  point  dipole  u  —  n  fi  (with  n  -  Q  (a)). 

0<  X  o(  0  1 

The  electrostatic  part  of  the  total  potential  is  the  sum  of  a  short 

o*.i  ir) 

ranged  part,  v  and  the  electrostatic  part  V  (2.10) 

We  have 

-  -,r  - 

V  (XaX  )-  Z0e  ^(XaX  ) 

-  Zo®  El!  <t>  ,  (0,*0  (3-2) 

where  4>  (0,^)  is  the  electrostatic  potential  generated  in  R0  by 

l 

particle  i  of  species  a. If  our  mixture  consists  only  of  charges  and  dipoles 
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E  in  Eq  (2.1)  verifies  E  --7  (with  V  -7  ).  Now 
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Using  the  relation 
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In  this  equation  Q—  Q(X-O) ,  where  Q(A)  is  the  partition  function 

Vs) 

defined  by  Eq.(2.5)  and  V  is  the  short  ranged  part  of  the  total  potential 
energy  (in  our  case  hard  core  interactions) . 

In  the  integral  in  (3.6)  the  first  term  is  zero  for  a  screening 
system,  in  which  the  field  and  potential  are  zero  at  the  system  boundaries. 
This  can  be  shown  using  Green's  theorem  and  Poisson's  equation. 
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Now,  since  the  media  particles  are  hard,  and  the  cavity  particle 
is  also  hard  and  of  arbritrary  shape, the  short  ranged  force  is 
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A  A 

where  a  (fl  R  ) 


r  (n  R  )  is  the  distance  of  closest  approach 
•> 

of  i  and  0,  which  depends  on  the  orientations  of  R  ,  and 


Using  (  3.3)  and  (3.9)  in  (3.8)  we  get 


where  <f>  is  the  electrostatic  potential  due  to  particle  a  (  3.3) 
*  v  t 

Separating  the  terms  in  which  a. -0.  .  we  get 
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The  self  term  is 
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where  the  sum  is  over  all  species  a  in  the  system  with  charge  z  e  and 
dipole  moment  /^-The  function  is  the  contact  probability  density 

between  the  particle  a  and  the  test  particle  O.It  is  a  function  of  the 

A 

angle  R  which  gives  the  position  of  a  relative  to  the  test  object  0. 

\j  A 

n(Rfl  )  is  the  unit  vector  along  the  center  to  center  distance  when  the 

A>  -A 

particles  are  in  contact. (R^  0o)  is  the  distance  of  closest  approach 
for  that  orientation  (see  figure  1) . 

The  mutual  term  is  more  complex 
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which  is  the  average  of  field  in  the  (n  (R0|))  direction  when  particle  1 
(of  species  a  )  is  in  contact  with  the  test  particle.  This  is  a 
complicated  term  because  it  involves,  a  3  body  correlation  function 
(particles  0,1,2)  .  However,  in  the  Onsager  limit,  in  which  the  fields 

outside  of  the  test  particle  are  completely  screened  out,  this  term  should  be 


be  zero, and  we  get 
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(3.14) 


where  <E0>S  *-s  gi-ven  hy  (3. 12).  In  the  case  of  a  spherical  test  object  in 


in  a  ionic  mixture, we  get 

t-\ 


<E  >  --4 n/(0  z  e) 

o  o 


(^r) 
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(3.15) 


IV -THE  ONSAGER  LIMIT 


In  his  classic  paper  of  1939,  Onsager14  considered  the 
problem  of  calculating  an  exact  lower  bound  for  the  interaction  energy 
V  of  a  neutral  system  composed  of  charged  hard  particles.  Onsager  solved 
the  problem  by  immersing  the  impenetrable  particles  in  an  infinite  conducting 
fluid,  thereby  obtaining  the  'Onsager  state'  for  which 


V  *  V 


(4.1) 


where  (the  Onsager  proper  energy  of  one  particle  of  species  a)  is  the 

self  energy  (relative  to  the  self  energy  of  the  original  charges)  of  the 

'Onsager  atom' .  This  object  is  the  original  hard  sphere  with  the 

original  charge  plus  that  surface  charge  induced  on  the  sphere  surface  by  the 

(£) 

conducting  boundary  condition.  Therefore  U^^is  the  self  energy  of  the 
(induced)  surface  charge. 

Onsager' s  approach  is  valid  also  for  point  charges  in  non-hard  core 
particles  (e.g.  the  one-component  plasma)  as  long  as  the  particles 
distribution  features  pair-excluded  regions  (i.e.  'effective'  hard  core)  10  19 


Recent  work  showed  10  13  20  25  that  the 'Onsager  state  ’provides  the 
exact  description  of  the  asymptotic  strong  coupling  limit  (e.g.  T-*0)  of  the 
solutions  of  both, the  mean  spherical  approximation  (MSA)  and  hypernetted  chain 


£  ft  \ 

(HNC)  integral  equations  for  arbritrary  hard  particles.  In  particular^  1 

which  is  the  sum  of  individual  terms,  one  for  each  particle  in  the  system, 
is  the  exact  strong  coupling  limiting  energy  as  obtained  from  the  solution 
of  the  MSA  and  HNC  equations . The  corresponding  asymptotic  direct  correlation 
functions  are  given  by  the  electrostatic  interaction  between  the  (induced) 
surface  charges  on  the  different  particles. 

Like  the  general  formalism  (Sec  II)  and  the  particular  MSA  formalism  (Sec 
V  below),  the  electrostatics  involved  in  the  'Onsager  State'  can  be 
developed  for  both  real  and  complex  charges. In  particular , relations  like 
(b.7)  in  appendix  B  follow  immediately  from  the  Onsager  limit  state 
due  to  the  simple  meaning  of  the  direct  correlation  function  as  stated  above. 

Consider  the  microfield  distribution  problem  in  the  Onsager  limit. 
Denoting  by  F(A)--kT  In  Q(A)  the  excess  configurational  free  energy  of 
the  system  when  the  imaginary  test  point  dipole  is  acting  on  the  test  object 
pointy  where  the  electric  microfield  is  to  be  considered,  then  Eq(2.4)  yields 


4 


i 


i 


a 


In  A(A)-  -0  [F(A) -F(A-O) ] 


(4.2) 


In  the  strong  coupling  limit,  the  system  approach  to  Onsager  state 
in  which  the  test  object  is  immersed  in  a  continuous  conducting  medium.  Then 
the  free  energy  difference  in  Eq  (4.2)  is 


In  A(A) | 


-0  [ U( A ) -U(A-O) ] 


(4.3) 


where  (see  Fig  2a)  U0(A)  is  the  self  energy  of  the  'Onsager  atom' 
corresponding  to  the  test  particle,  with  the  imaginary  dipole  iA//3  placed  at 
the  considered  point,  when  it  is  immersed  in  the  conducting  medium.  Since  we 
have  to  average  over  all  the  orientations  of  A , then  the  energy  difference 
(4.3)  is  proportional  to  the  square  of  the  dipole  strength  A: 

In  A(A)  |  -  A*7/9  <U  ^  (4.4) 

OtI  S  A  4  c  I*  /*■*' 


Here  U  ..  is  the  Onsager  proper  energy  (see  above)  of  a  hard  particle, 

having  the  shape  of  the  test  object  and  with  a  real  unit  dipole  which  is 

placed  at  the  test  point  (see  Fig.  2b).  <  denotes  the  average  over 

-n>* 

the  direction  of  the  dipole  m  in  a  frame  fixed  to  the  test  particle. 

Thus  the  Onsager  state,  as  well  as  the  strong  coupling  solutions  of  the 
MSA  and  HNC  equations,  predict  a  Gaussian  distribution  for  the  electric 
microfield  at  any  point  inside  the  test  object. The  second 
moment  depends  only  on  the  shape  of  the  particle  and  the  position  of  the 
test  point,  independent  of  the  charge  distribution  inside  the  test  object.lt  is 
given  by 

(4.5) 


m  -  <E  >  - ►  -  6kT  <U„  „  >  ^ 

u< 


r*'  -v 


The  Onsager  self  energy  for  a  dipole  m  at  the  center  of  a  sphere  of 

\/2)^/K  . 


radius  R  is 


For  example,  when  Che  unic  dipole  (i.e.  the  test  point  )  is  placed 
at  the  center  of  a  spherical  particle  of  radius  r  —a/ 2,  we  obtain 
regardless  of  the  original  charge  distribution  of  the  system, 

<E*  >  -3/(0  r5)  -2 4/(0  oJ) 

o  o>r«*<;*4 


(4.6) 


For  a  multicomponent  plasma  at  a  charged  point  of  charge  z  e 


the  Wigner-Seitz  ion  sphere  radius  is  a 


Vi 


,  where  <z>  is  the  mean 


electrovalence  .  Then  the  exact  second  moment  Eq  (3.1)  is  recovered: 


<E *”  >  -3/(0  &  )(<z>/zJ 


“4*  p  <z>/0zq 
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(4.7) 


Now  using  e/a^  as  the  unit  of  electric  field  (T-0 

<E*  >  -3/(T  )(  2 

0  0«uca<  0 
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(4.8) 


V  THE  MEAN  SPHERICAL  APPROXIMATION 


The  general  formalism  of  the  microfield  distribution  developed  in 
Section  III  is  suited  to  integral  equations  techniques.  The  appropiate 
Ornstein-Zernike  (O.Z.)  equation  for  the  fluid  in  the  presence  of  the  test 
particle  X  26 
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(5.1) 


where  c  (a,  £  -  0,1,.... s)  is  the  direct  correlation  function. 

In  Eq  (5.1)  we  have  included  the  test  object  as  member  of  a  new  species  with 
number  density  N0  /  V  ~Pa-  Eventually  we  will  take  the  limit,  No-»  1 
(or  in  the  thermodynamic  limit, p  —  0). 

The  Ornstein-Zernike  equation  can  be  closed  by  several  approximate 
forms  of  the  direct  correlation  function. The  simplest  of  such  closures  is 
the  mean  spherical  approximation  27,The  particles  are  hard  particles, 


so  that 


(5.3) 
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R><7 

‘1 


where  a  -  (a.  +a-v  )  /  2; a,  is  the  diameter  of  species  a 

<*•»,  *  r  * 

and  is  the  pair  potential. 

The  general  solution  of  Eq.(5.1)  with  the  closures  (5.2)  (5.3), 
for  arbritrary  multipolar  expansions  28  29 ,  shows  that,  because  of  the 
linearity  of  the  MSA,  the  R- integrals  in  Eq  (2.18)  are  proportional  to  the 
imaginary  dipole  strength. 


!  K  »*  V 
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(5.4) 


Then,  we  see  that  in  the  MSA  the  Fourier  transform  of  the  microfield 
distribution  is  Gaussian 

In  A(A)  -  -1/6  m^A  2 

(5.5) 


where  ■<  E  J  >,the  second  moment  of  the  distribution,  depends  on  the 
particular  multipolar  density  expansion  and  the  diameter  of  the  particles 
as  well  as  on  the  multipoles  and  diameter  of  the  test  object. 


We  can  take  advantage  of  the  availability  of  analytical  solutions 
of  the  MSA  for  some  particular  models  in  order  to  obtain  explicitly  their 
second  moments 

Lado®  has  recently  considered  the  microfield  problem  in  the  MSA 
for  a  one  component  plasma  of  point  charged  particles  in  a  neutralizing 
background  .  In  this  work  we  study  a  larger  class  of  systems: 

a)  Mixture  of  charged  hard  spheres  (primitive  model  of  electrolyte) 

b)  One  component  system  of  dipolar  hard  spheres  and 

c)  Mixture  of  charged  hard  spheres  and  dipolar  hard  spheres(non-primitive  model 
of  electrolyte) . 

For  these  systems  we  study  in  particular  the  strong  coupling 
limit  of  the  Onsager  states.  We  will  restrict  ourselves  to  the  special  case 
in  which  the  test  object  is  one  of  the  particles  of  the  system  and  will 
consider  only  the  microfield  at  the  sphere  center. 

a) PRIMITIVE  MODEL  OF  ELECTROLYTES 


The  system  is  a  mixture  of  s-1  species  of  hard  spheres  of  diameter 
a  with  a  point  charge  z^e  located  at  the  center  and  number  density 

P  (a  -  1 . s-1).  The  solvent  is  treated  as  a  structureless  continuum  with 

a  dielectric  constant  t • 

The  electric  field  felt  by  a  test  particle  of  type  7  with 


and  z^-zq  located  at  R  is, at  the  center  of  the  sphere 


(5.6) 


,  „  \  vi  t 

Hc=l 

The  correlation  function  which  involve  the  test  particle  (labelled  c>) 
has  the  form  (see  Eq  (2.17)) 


>  -»  <*\  lf)  N\  A 

f  (R  ;Q  )-  f°°°  (R)  +  f101  '(RJ  [iO  .RJ 


where  f  denotes  h  and  c  in  the  Ornstein  Zernike  equation  (5.1). 
Then  Eq  (2.7)  yields  (see  also  Eq  (2.18)) 


In  A(A)-4»/3(»/^^zidA'  dR  h“l<*'(R0l 
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As  shown  by  Lado9  for  plasmas ( see  Appendix  B),if  we  solve 
the  Ornstein  Zernike  equation  (5.1)  in  the  limit  ^-*0 
in  the  MSA  , then 


(5.7) 


(5.8) 


<E„,>  -v«  %•)  ■>^00 

(5.9) 
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Here  f°°°(R  )-f°°°  (R  )  is  the  correlation  function  for  the  unperturbed 

VH  v  ov'  **  0\'  r 

(without  the  imaginary  dipole)  fluid. 

Using  Eq  (5.9)  for  h101  and  taking  into  account  the 


electroneutrality  condition  we  obtain  Eq  (5.5)  with 


m  «-4»r/(/3  z  «) 

2.  o 


>  g  (*  J 

/  /  *  0*  o* 


(5.10) 


where  g  (a  )  is  the  contact  value  of  the  pair  distribution  function 
for  the  primitive  model  of  an  electrolyte  in  the  MSA. 

Comparison  of  this  second  moment  with  the  exact  result  given  by  (3.11) 
shows  that  the  MSA  has  only  the  self  contribution  <E*>  and  does  not  have  the 
mutual  term  <E*>  .This  is  to  be  expected,  since  the  MSA  is  an 

0  NN 

improved  Onsager  bound, should  not  have  this  term, which,  anyway 
vanishes  in  the  high  coupling  limit. 

If  we  use  the  MSA  result  of  the  contact  pair  correlation  function 
to  compute  the  limiting  form  of  m  in  the  high  coupling  regime, we  find 
that  it  diverges . This  is  a  consequence  of  the  well  known  inconsistency 
of  the  MSA, in  which  the  correlations  of  opositely  charged  ions  go  to 
-<*>  rather  than  to  zero. 

An  alternate  way  to  compute  the  Onsager  limit  of  the  second  moment  m^ 
is  to  include  the  test  particle  as  a  new  species.with  an  imaginary 
dipole. Our  system  consists  then  of  s-1  species  plus  a  new  kind  of 
particle  of  species  y  but  with  an  additional  imaginary  dipole. For 
this  system  the  energy  parameter 
aO 


J<r 


«o« 

dR  (R) 


is  known  in  closed  form  in  the  MSA29  34.The  result  is  proportional 
to  X  /p.  Taking  the  limit  p  -*  0,  the  energy  integrals  of  our 
eq.  (5.8)  converges. 

We  see,  then  ,  that  in  the  strong  coupling  limit  the  correct 
Onsager  limit  (4.6)  is  obtained. However  we  must  use  the  energy  route 
to  obtain  our  results,  because  of  the  well  known  inconsistency  of  the 


MSA. 


b)  HARD  DIPOLE  FLUID 


Consider  now  a  collection  of  hard  spheres  of  diameter  a 
with  an  embedded  point  dipole  /i  .We  compute  the  microfield  distribution 
at  the  center  of  these  particles . The  field  generated  by  a  particle  at 
the  center  of  the  test  particle  0  is 
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m 
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The  relevant  correlations  are 
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From  (2.7)  (and  also  from  (2.18))  we  have 
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As  shown  in  appendix  B  , the  radial  coefficients  hll0<^R  )  and  h112  (R  ,) 
rr  ’  os  o\  04  °> 

in  (5.12)  are  related  to  the  ordinary  (unperturbed)  radial  coefficients 
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(5.14) 


Again,  we  obtain  the  Gaussian  microfield  distribution  with  the  second 


(5.15) 


(5.16) 


(5.17) 

Indeed,  the  Onsager  limits  are  also  verified  in  this  case. 

c)  ION  DIPOLE  MIXTURE 

Consider  now  a  mixture  of  s-1  different  species  of  charged  hard 
spheres  of  diameter  ,  charge  z^e  and  number  density 
and  the  solvent  which  has  a  diameter  a  and  dipole  moment  n 

m  W 

The  field  at  the  center  of  the  test  particle  is  the  sum  of  the 
ionic  and  dipolar  contributions,  given  by  (5.6)  and  (5.11). 

On  the  other  hand  the  relevant  correlation  functions  now  are 

h*  (R  0  )  (5.7)  and  h  (X  ,X  )  (5. 12). The  microfield  is 

0«l  *  >  04  o  V 


moment  given  by 


m  -4b  /(30  a3  ) 

Z  2  4 


The  dipole-dipole  energy  parameter  is 


b  -*3»r  7 2/15  p  a3 
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and  which  can  be  calculated  solving  the  Wertheim32  equation 

(l+b,/3)2  -  (1-b,  /6)2  0  p,  M.2-d? 

(Ub‘/U)’  ^  ‘  ‘  1  * 


In  the  strong  coupling  limit  d2  -»  ®  and  b^  -*6 . 


If  the  test  particle  is  an  ion  then  we  again  use  the  method 


of  defining  an  extra  species  A  of  diameter  o, charge  z  e  and 
imaginary  point  dipole  of  strength  -i\/P  .  We  now  let  the  density 
of  this  hypothetical  component  go  to  zero.  Eq .  (5.18)  yields 
In  A(A)  —  4/(3/J)  IdA'  [b^/ZfZ  +b  (A) /2a3  ] 

Jo  (5.21) 

The  calculation  of  b  and  b  as  functions  of  A  is  identical 

»  z. 

to  the  case  of  the  real  dipole, and  is  described  elsewhere33  34.The 
important  point  is  that  in  the  MSA,  b(  (A)  and  b^(A)  are  linear 
functions  of  A  .Therefore,  we  obtain  again  the  Gaussian  distribution  (5.5). 


VI -CONCLUSIONS 


In  this  paper  we  discuss  the  microfield  of  a  ionic  solution  in  a 
molecular  solvent. We  extend  the  Morita-Iglesias  formalism 
to  include  angular  dependent  interactions . 

The  mean  spherical  approximation  yields  a  simple  answer  for  the  I 

microfield  in  all  the  cases  that  were  investigated: The  primitive  model 
of  ionic  solutions,  the  pure  dipolar  fluid  and  the  ion-dipole  mixture. 

For  all  of  these  cases  the  distribution  is  a  Gaussian. The  width  of  1 

these  Gaussians  (which  is  the  inverse  of  the  second  moment)  can  be  obtained 
in  two  forms : 

i)  in  terms  of  the  contact  probability  .  This  form  agrees,  in  the 
Onsager  limit,  with  the  limiting  form  of  the  exact  value  .also  derived 
in  this  paper. 

ii)  in  terms  of  the  excess  energy  parameter.  I 

An  interesting  situation  arises  here:  In  the  Onsager  limit,  the  MSA 

contact  probability  of  the  ions  diverges.  So  that  although  formally  the 
MSA  gives  the  correct  limiting  expression  of  the  second  moment  in  terms 
of  the  contact  probabilities,  the  limit  itself  is  incorrect. 

The  energy  route  (ii),  however,  leads  to  the  correct  limiting  value  of  the 
second  moment  in  the  Onsager  limit. This  is  probably  a  consequence  of  the 
well  known  inconsistency  of  the  MSA. 

Sifmilarly  to  the  situation  for  plasmas9  the  MSA  should  serve 
as  a  starting  point  for  a  more  accurate  integral  equation  treatment. In 
this  treatment  the  Onsager  gaussian  result  (4. 4), (4. 5)  is  the  exact 
high  coupling  limit. 


APPENDIX  A 


ONSAGER  ENERGY  FOR  A  DIPOLE  IN  A  SPHERICAL  CAVITY 


Following  Rosenfeld  and  Blum12  13  we  write  the  Onsager 
energy  of  an  arbitrarily  charged  object  immersed  in  a  conducting  medium 


(a .  1) 

--[self  energy  of  the  surface  charge  distribution  o(s)] 

Here  ^  (s)  is  the  potential  due  to  the  induced  charges  at  the 

c>t 

surface  E.For  a  spherical  object  with  a  unit  dipole  at  the  center  (and 
pointing  in  the  direction  of  the  z  axis)  we  have 


and 
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where  9  is  the  angle  between  the  normal  to  the  surface  and  the  z  axis. Then 
is  given  by 


U-  - 
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Finally,  from  (a. 4)  and  (4.5)  we  get  (4.6). 


APPENDIX  B 


CORRELATION  FUNCTIONS  IN  THE  MSA 
Consider  the  Ornstein-Zernike  equation  (5.1)  in  the  limit  po  -0 
We  use  the  invariant  expansion  (2.17)  for  all  the  correlation  functions. 
Then, using  the  transforms28  28 
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where  j^(kR)  is  the  spherical  Bessel  function  of  order  l.and  F  is  a  generic 
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notation  for  either  C  or  H  (transforms  of  c  or  h).The  Ornstein-Zernike 
equation  is 
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(b .  2) 


The  Ornstein  Zernike  equation  can  be  transformed  to  real  space :Using 
the  inverse  transform  of  (b.l) 


oa 


a ) ±\ki 

^f>;X  ^  K  J  ~(vr)wt  iZ  \  /< 


rZ) 

eos{«>4j 


I  I/I 


(b .  3) 


where  the  upper  sin  kR  is  for  m+n  even  ,and  the  lower  cos  kR  is  for  m+n 
odd. 

We  analize  now  the  case  of  the  microfield  at  the  center  of  one 


of  the  ions  of  the  mixture,  and  at  the  center  of  a  hard  dipole .which 
could  be  in  a  ion-dipole  mixture  or  a  pure  dipolar  fluid, 
a) TEST  ION  IN  THE  PRIMITIVE  MODEL 

Using  (5 . 7) , together  with  (b.l)  and  (b.3)  in  the  Ornstein-Zernike 


equation  (b.2)  yields 
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The  observation  is  that  the  MSA  closure  relations  (5.2)  and  (5.3) 


yield 
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Then,  defining 
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for  all  R,  and  using  the  fact  that  roo^  {  U  ) 
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so  that  if  0  Cft\  verifies  (b.4),  then  J  l  defined  by  (b.9) 

ad  i  o  d  ’ 

should  verify  (b .  5)  .Therefore ,  f^l(^R  )  as  given  by  (b.8)  is  the  solution 

of  the  problem. We  also  notice  that  f^° ^ftyis  equal  to  f^0(<jfor  the 

bulk  unperturbed  system  (A-0)  when  the  test  particle  belongs  to  species  y. 


(b.9) 


5  4 


b)TEST  DIPOLE  IN  AN  ION- DIPOLE  MIXTURE 

For  this  case  the  one  dimensional  projection  of  the  Ornstein-Zernike 

equation  yields  (y-r.A)  {+>4 
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(b.12) 


The  MSA  closure  conditions  are 


r-J(  .1 e.?n, 


(b. 13) 


(b.14) 


/5  /6y  __  J  cJ  if  ^  4  /  toii;/  '  ^  +  <§5i  //)  J 

b  ~  ( 

with  for  7—r  and  for  7-A.But  we  observe 

that  for  y-r  equations  (b. 10-14)  are  the  MSA  equations  for  for  the 
unperturbed  system  (A-O).This  problem  is  solved  elsewhere29  34 
When  7-A  we  have 


f <T~ '**"*/  fJ 

CTo*  (  (*«.+)  ■*  (i  °J;  (<SJ  J 

^MS  V 


(b. 15) 


The  correlation  functions  have  coefficients 

)-  -(V/3  b^1  (R  )  s-1 

h“of>U  )-  -(V/3  /*)  h^°  (R  ) 
h^2W(R  )-  -(A/0  M/>  (R  ) 

(b . 16) 

The  case  of  the  pure  dipolar  fluid  is  obtained  from  this  result  by 
setting  the  ionic  concentration  equal  to  zero. 
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FIGURE  CAPTIONS 

Figure  1: 

Distance  of  closest  approach  between  a  convex  test  particle  and  hard 
spherical  molecules. 

Figure  2: 

Self  energy  of  Onsager  objects  immersed  in  a  perfect  conductor. 


